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AN  ATTITUDE  DETERMINING  SYSTEM  FOR  SPINNING  ROCKETS 

ABSTRACT 

An  attitude  sensor  for  spinning  rockets  and  the  associated  digital 
computing  procedure  have  been  developed  and  tested.  The  measuring 
system  consists  of  a  commercially  available  flux  gate  magnetometer  and 
a  single  solar  cell  with  a  fan-shaped  acceptance  beam.  The  computing 
procedure  is  designed  to  improve  approximations  of  the  rocket  aspect 
angles  by  successive  differential  corrections  which  are  derived  by  a 
least  squares  treatment  of  an  overdetermined  system  of  condition  equations. 

Both  the  measuring  system  and  the  computing  methods  have  been  found 
to  provide  a  practical  approach  for  determining  rocket  aspect.  The  system 
is  estimated  to  have  been  accurate  to  between  one  and  two  degrees  in  its 
initial,  test  and  to  have  a  potential  accuracy  of  better  than  one  degree. 
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the  Jacobian  matrix  ( - g — - - ) 

vector  in  the  direction  of  the  axis  of  the  magnetometer 
the  numerator  of  Equation  (l8) 
a  rotation  matrix 

vector  from  the  origin  toward  the  sun 
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Subscripts 


computed  value 


end  of  reduction  interval 
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observed  value 


refers  to  vector  C 
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Superscripts 
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INTRODUCTION 


The  purpose  of  the  efforts  reported  here  was  to  obtain  moderately 

accurate  (±  1°)  measurements  of  sounding  rocket  attitude  with  a  minimum 

of  construction  expense,  time,  and  calibration.  Because  of  the  ready 

availability  of  a  reliable  and  easily  used  flux  gate  magnetometer,  a 

large  proportion  of  sounding  rockets  incorporate  such  an  instrument  for 

qualitative  aspect  information,  i.e .,  orientation  with  respect  to  the 

1  2* 

geomagnetic  field  vector  only  9  .  Israel,  et.  al.,  have  reported  the 

3 

determination  cf  complete  aspect  information  from  magnetometers  only  . 

Solar  aspect  devices  have  also  been  used  for  qualitative  aspect  measure  - 

1  5 

raents  with  respect  to  the  rocket-sun  vector  .  Schemes  involving  one 
or  more  magnetometers  and  calibrated  or  uncalibrated  photo -detectors 
have  been  reported  for  acquisition  of  complete  orientation  information  9  . 
The  system  described  herein  has  the  advantage  that  only  an  uncalibrated 
solar  cell  or  photomultiplier  need  be  added  to  the  normally  present 
aspect  magnetometer  in  order  to  obtain  complete  attitude  information  for 
a  spinning  rocket.  The  magnetometer  is  factory  calibrated  and  yields 
the  angle  between  the  rocket  axis  and  the  geomagnetic  field.  The  solar 
cell  is  collimated  to  a  narrow,  fan -shaped  beam,  and  the  phase  of  the 
solar  cell  pulse  relative  to  the  magnetometer  signal  is  used  to  complete 
the  attitude  determination. 

A  digital  computing  method  has  been  developed  for  the  reduction  of 
the  flight  measurements  to  rocket  aspect.  This  problem  requires  the 
evaluation  of  a  set  of  three  parameters  to  completely  specify  the 
orientation  of  a  rigid  body.  A  convenient  set  of  parameters • consists 

Q 

of  the  Eulerian  angles  ,  two  of  which  determine  the  attitude  of  the 
body  while  the  third  describes  the  rotation  of  the  body  about  its  symmetry 
axis.  By  combining  certain  simplifying  assumptions  with  the  observations 
recorded  for  the  solar  cell,  it  becomes  possible  to  express  the  third 
parameter  as  a  function  of  the  other  two  Eulerian  angles.  Hence,  the 
problem  reduces  to  the  determination  of  two  angles  which  vary  with  time. 

*  Superscript  numbers  denote  references  which  may  be  found  on  page  47. 


Moreover,  these  angles  completely  describe  the  aspect  of  the  rocket. 
Since  the  measured  quantities  are  nonlinear  functions  of  the  Eulerian 
angles,  it  is  convenient  to  evaluate  these  angles  by  a  fitting  procedure 
in  which  a  set  of  approximations  is  improved  by  a  series  of  differential 
corrections.  The  latter  are  obtained  from  a  least  squares  treatment  of 
an  overdetermined  system  of  condition  equations  which  arc  derived  from 
the  measurements  of  the  aspect  system. 

The  method  of  solution  was  tested  on  actual  field  data  for  a  Nike- 
Apache  missile  launched  at  Fort  Churchill,  Canada.  Both  the  measuring 
system  and  the  computing  procedures  were  found  to  provide  a  practical 
approach  for  determining  rocket  aspect.  Generally,  the  computations 
converged  to  the  correct  solution  on  the  BRLESC  computer  in  approximately 
0.6  seconds  per  point.  It  Is  estimated  that  the  aspect  system  provided 
an  accuracy  of  better  than  two  degrees  in  its  initial  test  and  further, 
that  its  potential  accuracy  is  one  degree  or  better. 

INSTRUMENTATION 

Figure  1  is  a  perspective  diagram  of  the  magnetometer  sensor  and 
solar  cell  within  the  rocket  payload.  The  aspect  magnetometer  was  a 
flux  gate  instrument  capable  of  about  1  percent  accuracy.  It  was  glued 
to  a  fiberglass  deck  plate  after  being  located,  such  that  it  was  normal 
to  the  plane  containing  the  optical  aperture  and  the  symmetry  axis  of 
the  rocket.  The  output  of  the  magnetometer  in  volts  (E)  is  proportional 
to  the  component  of  the  geomagnetic  field  parallel  to  the  sensor  axis. 
That'  is, 

E  -  4.0  |h|  cos  X  ,  (1) 

where  X  is  the  angle  between  the  geomagnetic  field  vector  H  (gauss) 

— * 

and  the  sensor  positive  axis  M.  The  solar  cell  was  a  commercial  unit 
5  mm  x  5  mm,  mounted  one  inch  from  a  longitudinal  aperture  In  the  fiber¬ 
glass  skin  of  the  rocket.  The  aperture  was  covered  by  a  fused  quartz 
window  to  pi  otect  the  solar  cell  from  aerodynamic  effects,  but  this  was 
probably  not  necessary. 
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For  a  rocket  with  spin  frequency  much  greater  than  pitch  frequency, 
the  transverse  aspect  magnetometer  produces  a  modulated  sinusoidal  out¬ 
put  signal,  the  instantaneous  amplitude  of  which  defines  a  cone  about 
the  geomagnetic  field-  As  shown  in  Figure  2,  one  of  the  generating 
elements  of  this  cone  coincides  with  the  rocket  longitudinal  axis.  In 
the  aspect  system  described  here  a  solar  cell  with  a  fan -shaped  acceptance 
beam,  determined  by  the  slit  in  the  rocket  skin,  produces  a  pulse  as  it 
sweeps  past  the  sun.  Then  the  phase  of  this  sun  pulse  relative  to  the 
magnetometer  signal  is  a  function  of  the  position  of  the  rocket  on  the 
previously  defined  cone.  Figure  2  illustrates  the  relationships  between 
the  geomagnetic  field,  the  cone  defined  by  the  aspect  magnetometer,  and 
the  solar  cell  acceptance  beam. 

The  solar  cell  and  the  slit  were  both  5  *nm  vide,  and  the  cell  was 
located  25  mm  behind  the  slit.  Thus,  the  angular  aperture  of  the  solar 
cell  in  the  horizontal  plane  was  23°,  but  the  cell  was  fully  illuminated  onl, 
instantaneously ,  at  the  center  of  the  pulse.  This  produced  an  approxi¬ 
mately  triangular  waveform,  later  differentiated  and  bandwidth  limited. 

The  length  of  the  3iit  was  3  inches,  yielding  an  angular  coverage 
of  +  66°  to  -  i*l°,  since  the  solar  ceil  was  not  centered  vertically. 

The  ooiid  angle  coverage  was  3*1  n  steradians,  adequate  for  the  present 
application.  This  could  be  readily  extended  to  3*8  n  steradians  with 
the  sume  slit,  and  two  solar  cells  in  parallel,  located  behind  the  upper 
and  lower  edges  of  the  3"  aperture. 

The  magnetometer  signal  was  direct  coupled  into  a  standard  IRIG 
10.5  kc  cubcarrier  oscillator  (SCO)  of  an  FW/FM  telemetry  system,  as  shown 
in  the  electrical  schematic  (Figure  3)*  The  common  base  amplifier  used 
with  the  solar  cell  produces  approximately  5  volts  signal  for  solar 
illumination  of  normal  incidence  above  the  atmosphere.  Since  the  pay- 
load  construction  was  nearly  complete  when  the  solar  cell  was  added, 
the  sun  signal  was  merely  AC  coupled  into  the  same  SCO  used  for  the 
magnetometer.  The  delay  of  the  sun  signal  relative  to  the  magnetometer  . 
signal,  at  the  discriminator  output,  was  approximately  0.3  milliseconds 
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FIGURE  3  -  ELECTRICAL  SCHEMATIC 


for  spin  rates  of  5  to  10  rps,  and  was  considered  negligible  for  this 
application.  For  automatic  data  reduction  separate  channels  for  the 
two  functions  would  be  desirable,  with  appropriate  correction  for  the 
greater  differential,  delay. 

Two  typical  oscillograms  from  Nike -Apache  rocket  No.  14.36  are 
shown  in  Figure  4.  The  slowly  varying  sinusoidal  waveform  is  the 
magnetometer  output  signal,  and  the  transient  sun  pulse  is  superimposed 
on  it.  The  differentiating  effect  of  the  AC  coupling  and  the  bandwidth 
limiting  of  the  discriminator  filter  produced  the  sun  pulse  shapes 
illustrated.  The  magnetometer  signal  is  essentially  undistorted.  Sample 
A  is  from  early  in  flight  while  the  rocket  was  nearly  vertical.  The 
phase  of  the  sun  pulse  relative  to  the  magnetometer  signal  remained 
nearly  constant  during  this  interval.  Sample  B  occurred  during  rocket 
turnover  upon  re-entry  into  the  effective  atmosphere.  The  sun  pulse 
can  be  seen  to  shift  about  l60°  in  phase  during  2.6  seconds  of  flight 
time.  The  sun  was  outside  the  solar  cell  field  of  view  during  the  first 
three  and  the  last  two  of  the  cycles  shown.  This  particular  flight 
produced  poor  quality  telemetry  data.  There  was  complete  loss  of  sig¬ 
nal  from  37  seconds  to  265  seconds  and  frequent  short  dropouts  occurred 
from  265  seconds  to  rocket  impact.  Even  so,  the  aspect  system  produced 
usable  data  at  all  times  for  which  the  sun  signal  could  be  detected. 

DEVELOPMENT  OF  THE  COMPUTING  PROCEDURE 

The  method  of  solution  used  in  reducing  the  recorded  data  to  missile 
aspect  is  predicated  upon  the  assumption  of  a  negligible  change  in 
rocket  attitude  over  the  time  interval  of  a  single  revolution  of  the 
missile.  In  fact,  for  rapidly  spinning  rockets,  we  may  assume  that 
the  aspect  remains  constant  throughout  each  revolution,  but  varies  from 
one  revolution  to  another.  While  this  assumption  is,  perhaps,  less  valid 
for  slowly  spinning  missiles,  it  is  still  a  reasonable  approximation  for 
at  least  the  free  flight  portion  of  the  trajectory.  Hence,  we  consider 
successive  pulses  from  the  solar  cell  to  exactly  define  the  time  interval 
of  one  complete  revolution  of  the  rocket  and  assume  its  aspect  to  remain 

15 


FIGURE  4  “  SAMPLE  OSCILLOGRAMS  FROM  NIKE  -  APACHE  14.36  —  ASPECT 

MAGNETOMETER  OUTPUT  VOLTAGE  VS.  ROCKET  FLIGHT  TIME 


constant  between  pulses.  In  addition  to  the  solar  pulses,  only  two 
observations  from  the  magnetometer  are  required  to  compute  a  pair  of 
angles  that  completely  describe  the  attitude  of  the  missile.  Normally, 
ten  to  fifteen  readings  are  made  between  successive  pulses  to  provide 
considerable  overdetermination  for  the  solution.  The  method  of  computation 
consists  of  a  fitting  procedure  in  which  a  set  of  approximations  for 
the  aspect  angles  is  improved  by  successive  differential  corrections. 

These  corrections  are  obtained  from  a  least  squares  treatment  of  the 
overdetermined  system  of  condition  equations  that  are  derived  from  a 
Taylor  expansion  about  a  point  consisting  of  approximations  for  the 
aspect  angles.  All  second  and  higher  order  terms  of  the  expansion 
have  been  neglected. 

The  solution  will  be  described  for  a  single  revolution  of  the  rocket. 
The  computing  method  presented  here  is  then  merely  repeated  for  each 
revolution  of  the  missile  over  its  entire  trajectory.  We  define  t  and 
t^  to  be  the  times  of  successive  solar  pulses  and  hence ,  based  on  our 
assumptions,  the  initial  and  final  times  of  one  complete  revolution  of 
the  rocket.  The  symbol  t  will  denote  the  times  of  all  observations 
which  are  recorded  for  the  magnetometer  within  the  interval  [tQ,t^]. 

The  subscript  i  ranges  from  1  to  k  where  k  is  the  number  of  observations 
in  the  interval. 

In  the  development  of  this  solution,  it  is  convenient  to  employ 
two  rectangular  coordinate  systems,  one  fixed  with  respect  to  the  earth 
and  the  other  fixed  with  respect  to  the  orientation  of  the  missile  body. 
Referring  to  Figure  5,  the  two  systems  are  defined  as  follows. 

1.  The  XYZ -coordinate  system  is  a  right-hand  system  with  the 
origin  on  the  earthTs  surface  at  the  launcher.  The  X-axis  is  positive 
south,  the  Y-axis  is  positive  east,  and  the  Z-axis  is  normal  to  the, earth's 
surface  at  the  launching  site.  This  will  be  referred  to  as  the  space 
coordinate  system. 


2C  The  X'Y'Z*  -  coordinate  system  is  a  right-hand  system  with 
the  same  origin  as  the  XYZ-system.  The  2 1 -axis  is  parallel  to  the  symmetry 
axis  of  the  rocket  body  while  the  Xf  -axis  is  positive  in  the  direction 
of  the  magnetometer  axis.  The  Yf -axis  is  chosen  sc  as  to  complete  a 
right-hand  system.  This  will  he  referred  to  as  the  body  coordinate 
system. 


In  our  instrumentation  system,  the  YT  -axis  passes  through  the  slit  which 
has  been  cut  in  the  skin  of  the  rocket  for  the  solar  cell.  This  is  a 
simplifying,  but  not  a  necessary,  condition  for  our  method  of  solution. 

It  is  convenient  to  specify  the  orientation  of  the  body  coordinate  system 
with  respect  to  the  space  system  by  means  of  the  Eulerian  angles  0,  i|r, 
and  0  as  indicated  in  Figure  5*  These  angles  specify  the  direction  and 
magnitude  of  three  rotations  which,  when  performed  in  a  specific  sequence, 
will  result  in  a  transformation  from  one  of  the  rectangular  coordinate 
systems  to  the  other.  The  angles  0,  0,  and  i|r  are  generated  by  rotations 
about  the  2-axis,  the  line  of  nodes  OP,  and  2’ -axis  respectively.  All 
are  positive  for  a  right-hand  rotation.  Let  the  notation  R^  (7)  indicate 
the  matrix  performing  a  rotation  through  an  angle  7  and  the  axis  of 
the  frame  of  reference,  such  that  the  angle  is  positive  for  a  rotation 
in  the  right-hand  direction.  The  symbol  j  is  equal  to  1,  2,  or  3  according 
to  whether  the  rotation  is  about  the  X,  Y,  or '2  axis,  respectively.  In 
this  notation,  a  transformation  from  the  space  system  to  the  body  system 
can  be  accomplished  by  three  rotations  as  follows:  l)  R^(0),  2)  R^(9), 
and  3)  To  transform  from  the  body  system  to  the  space  system 

requires  the  following  sequence  of  rotations:  l)  R^-^);  2)  R^(-0),  and. 

3)  R3(-4). 


On  the  basis  of  previous  assumptions,  0  and  0  are  considered  to  be 
constant  over  the  time  interval  between  successive  solar  pulses,  but 
are  assumed  to  vary  from  time  interval  to  time  interval.  That  is  to  say, 
0  and  0  are  both  approximated  by  step  functions  which  change  values  on 
the  discrete  set  of  discontinuities  defined  by  the  solar  cell  pulses. 

The  solution  is  designed  to  evaluate  these  approximating  functions  step 


\ 


by  step.  Rocket  aspect  will  be  determined  when  0  and  0  are  known;  but 
if  desired,  the  results  may  easily  be  expressed  in  terms  of  the  more 
familiar  azimuth  and  elevation  angles.  In  this  method  of  solution,  there 
exist  ambiguities  which  can  normally  be  resolved  without  difficulty.  We 
introduce  constraints  to  force  convergence  to  the  correct  choice  of  four 
possible  solutions.  To  accomplish  this,  the  quadrant  of  the  angle  0  and 
the  direction  of  the  spin  of  the  missile  must  be  known.  The  latter  is 
usually  predetermined,  but  can  be  measured  if  necessary,  while  the  former 
is  readily  available  through  a  proper  interpretation  of  the  trend  in  the 
amplitude  of  the  magnetometer  data.  The  necessary  constraints  follow: 

1)  0  ^  0  <  tt/2  when  the  missile  is  nose -up, 

2)  rr/2  <  0  ■<:  tt  when  the  missile  is  nose -down, 

■ 

3)  V  >  0  for  right-hand  rotation, 

k)  \|r  <  0  for  left-hand  rotation, 

where  is  the  spin  rate  of  the  rocket.  If  the  results  are  desired  in 
terms  of  the  more  familiar  angles  of  azimuth  and  elevation,  we  note 
from  Figure  6  that 


e  »  tt/2  -  0  , 

(2) 

(  0  -  tt/2  when  rr/2  '<  0  <:  2n  , 

(3) 

/  0  +  3rr/2  when  0  ^  0  ^  tt/2 

where  a  is  defined  as  the  azimutn  measured  counterclockwise  from  the 
X-axis  and  restricted  to  a  range  of  from  zero  to  2rr  radians,  while  e  is 
the  elevation  angle  with  respect  to  the  XY -plane. 

In  addition  to  our  previous  assumptions,  we  consider  (|r  to  be  constant 
between  successive  solar  cell  pulses  so  that 

'  ,  2tt  /.v 

*  ■ 1  nf  --v  >  ■  w 

If  is  the  value  of  the  Eulerian  angle  |  at  time  t^  it  follows  that 
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(5) 


+i  =  +o  +  *(ti  "  to)  ' 


where  ijr^  is  the  unknown  value  of  t|r  at  time  tQ.  Shortly,  it  will  be 
shown  that  ijfQ  may  be  expressed  as  a  function  of  8  and  0.  Hence,  the 
probl  m  reduces  to  one  of  determining  the  step  functions  which  approximate 
0  and  0. 


In  our  development,  we  require  notation  for  the  direction  cosines 
of  several  vectors  in  both  coordinate  systems.  These  will  be  denoted 
by  (l,  m,  n)  in  the  space  coordinate  system  and  (l1,  m',  nT)  in  the 
body  system.  In  addition,  the  direction  cosines  will  be  further  identified 
by  the  subscripts  M,  C,  S,  and  H  which  correspond  to  the  following  vectors: 

l)  M  is  defined  as  a  vector  in  the  direction  of  the  axis  of 
the  magnetometer. 

2;  C  is  a  vector  lying  in  the  X!YT -plane  and  directed  from  the 
origin  of  the  coordinate  system  through  the  slit  cut  in  the 
missile  skin  for  the  solar  cell. 

3)  £  is  a  vector  from  the  origin  of  the  coordinate  system 

toward  the  sun. 

H  is  the  earth1 s  magnetic  field  vector. 

In  the  simplest  approach  to  the  problem,  the  vectors  S  and  H  may  be 
considered  constant  in  direction  over  the  entire  trajectory  of  the  rocket. 
For  greater  accuracy,  however,  S  can  be  varied  with  time  and  H  with 
geometrical  position  provided  that  missile  position  is  known  as  a  function 
of  time- 


In  the  development  of  the  solution,  the  direction  cosines  for  the 

—f 

vectors  II  and  S  will  be  considered  as  known  in  the  space  reference  frame. 

— > 

The  direction  of  H  may  be  determined  either  by  direct  measurement  in  the 

launching  area  or  from  series  expansions  using  published  spherical 

9 

harmonic  coefficients  while  the  direction  cosines  of  S  may  be  readily 
computed  from  ephemeris  of  the  sun^.  Hence,  we  regard  ‘(l^,  m^,  n^)  and 
(l  ,  mc,  nc)  as  known.  The  direction  cosines  of  the  vectors  M  and  C  are 

ODD 

constant  and  known  in  the  body  coordinate  system,  but  vary  with  time  in 


22 


the  space  system*  We  may  convert  from  the  body  to  the  space  system  by 
applying  a  transformation  matrix  which  is  the  triple  product  of  three 


separate  rotations.  Let  this  time  variable  matrix  be  detnoted  by 
Then, 


where 


It  follows  that 


and 


A1  »  R3(^)Rl(-0)R3(^.)  , 

(cos  0  -sin  0  O'' 
sin  0  cos  0  0 
0  0  1/ 

(10  0 

0  cos  0  -  sin  6 
0  sin  0  cos 

/cos  i|r  -  sin  O' 

|  sin  fj,  cos  ^ 

3  0 


M 

“Mi 

P*  I 


Ai 


“M 

*k 


/ 1  \ 

/ 1*  \ 

Ci  1 

1  C 

m„ . 

=  A, 

mi 

Ci 

i 

C 

J 

In'  / 

(6) 


(7) 


(8) 


(9) 


(10) 


(H) 


The  subscript  i  for  the  direction  cosines  in  the  space  coordinate  system 
simply  indicates  that  these  quantities  vary  with  time,  and  therefore, 
refers  to  the  particular  time  t^. 
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By  virtuve  of  the  definitions  of  the  body  coordinate  system  and 


the  vectors  M  and  C  we  have 


cos  ce 


=  sin  oj 


n  •  /  \  0 


where  cjo  is  defined  as  the  angle  between  C  and  X'-axls.  For  our  particular 

— >  s 

Instrumentation  system,  C  corresponds  to  the  Yf-axis  so  that  co  *»  rt/2  and 


"c  i"  x 


Y'cj  \°! 

With  this  simplification,  it  immediately  follows  that  the  direction 
cosines  (1^,  ^ii^  are  ®^ven  ^  the  elements  of  the  first  column 

of  the  matrix  while  (l  ,  n^)  ai'c  obtained  from  Hie  second 

column  of  the  same  matrix.  Hence.,  we  have 

1^  «  cos  0  cos  4'^  -  sin  0  cos  0  sin  ^  , 

«=  sin  0  cos  ^ .  +  cos  0  cos  6  sin  ,  (J 


nk,.  »  sin  0  sin  , 
Mi  i 
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1^  -  -  (cos  0  sin  ^  +  sin  0  cos  0  cos  , 

m,,,  =  cos  0  cos  0  cos  -  sin  0  sin  1/.  ,  (16) 

Li  l  j» 

m-,.  -  sin  0  cos  i.  . 

Ci  i 


It  will  now  be  shown  that  fQ,  the  value  of  at  time  t^,  may  be 
expressed  in  terms  of  the  unknowns  0  and  0.  The  vector  C  is  in  the 
direction  of  the  Y* -axis .  Referring  to  Figure  2,  it  is  apparent  that 

—4 

the  vector  3  must  lie  in  the  Y'Z* -plane  at  t  ,  a  time  at  which  a  pulse 
is  recorded  for  the  solar  cell.  Moreover,  the  vector  M  is  in  the 

—4 

direction  of  the  X' -axis  so  that  M  and  3  are  orthogonal  at  tQ.  The 
condition  for  orthogonality  is 

Vw>  +  “sV  '  V*>  ‘  0  •  <17> 


Substituting  from  Equations  (it)  with  i  =  0  and  solving  for  yields 


tan 


-1 


lg  cos  0  +  mg  sin  0 

1~  cos  0  ^irT0  -  rig  cos  0  cos  0  -"n^  sin  8 


(18) 


If  we  restrict  to  lie  between  0  and  2tt  radians,  there  are  two  possible 


■iOlUolUUb  lUi 


itch  satisfy  Equation  (18).  However,  we  have  an  added 
condition  which  permits  the  complete  resolution  of  this  ambiguity.  Since 
the  solar  cell  is  activated  at  time  tQ,  it  is  apparent  that  the  sun  is 
then  contained  in  the  solar  cell  acceptance  beam,  i.e.,  in  the  Y 1 Z ' -plane 
and  on  the  positive  side  of  the  X»Z' -plane  (Figure  l).  Hence,  the  angle 

>  I  "4  l 

between  S  and  C  must  be  equal  to  or  less  than  n/2  radians  so  that 

y-co +  Vco +  Vco $  0  ■  (19) 

Combining  Equations  (l6)  with  the  Inequality  (19)  yields 


sin  i|f  (-  lq  cos  0  -  nig  sin  0) 

+  cos  ^Q(  -  lg  sin  0  cos  9  +  nig  cos  0  cos  0  +  ng  sin  6)  S7  0, 


(20) 
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If  we  define  N  and  D  to  be  equal  respectively  to  the  numerator  and 
denominator  of  the  fraction  in  Equation  (18),  Inequality  (20 )  may  be 
written 

-  (N  sin  +  D  cos  tQ)  ^  0  .  (21 ) 

Combining  Equation  (l8)  with  Inequality  (21)  yields 

-  D  sec  *Q  »  0  ,  (22) 

-  N  esc  f  ^  0  .  (23) 

Inequalities  (22)  and  (23)  uniquely  determine  in  which  quadrant  i|»o  lies. 
The  results  are  tabulated  in  Table  I. 


Range  of  N 
N  <C  0 
N  <  0 
0  <  N 
0  <  K 


TABLE  I 

DETERMINATION  OF  THE  RANGE  FOR  f  IN  RADIANS 


Range  of  D 
D  <  0 
0  ■<  D 
0  <  D 
’£>'<  Q 


Range  of 
0  ^  t  <  n/2 
rr/2  "<:  i  <  tt 
tt  ^  tQ  <  3tt/2 

>tt/2  *Sf  *  *r-  />rr 


It  is  worth  noting  that  can  sti1 1  be  expressed  as  a  function  of 
0  and  0  when  to  is  not  equal  to  n/2  radians.  The  procedure  is  completely 
analogous  to  the  above  derivation.  The  results  will  be  presented  without 
development.  We  have  that 


,  j  N  sin  to  +  D  cos  to 

tan  v  *  — — - — - _  _  .  ) 

To  D  sin  to  -  N  cos  to 

with  the  conditions  that 

J  -  sec  +0(D  sin  to  -  N  cos  to)  ^  0  , 

-  esc  +  (D  cos  to  +  N  sin  to)  ?  0  . 
o 


(2k) 

(25) 

(26) 
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As  before,  these  permit  a  unique  determination  for  i|r  in  the  range  from 
0  to  2n  radians. 

We  now  consider  the  magnetometer  measurements  from  which  we  may 
derive  an  overdetermined  system  of  condition  equations.  These  measurements 
are  proportional  to  the  component  of  the  geomagnetic  field  vector  in  the 
direction  of  the  magnetometer  axis,  that  is,  the  X1-  axis  of  the  body 

coordinate  system.  After  normalization,  the  data  constitute  measurements 

— ♦  — * * 

of  cos  X^  where  X^  is  defined  to  he  the  angle  between  the  vectors  M  and  H 
at  time  t^.  We  may  then  write  the  system  of  equations 

¥mi  +  Vta  +  Vmi  ’  003  h  >  (27) 

for  all  values  of  i.  In  this  system,  the  direction  cosines  (l^,  m^,  n^) 
are  considered  known.  Referring  to  Equations  (5)^  (15)  a^d  (l8)  we  note 
that  (1^,  n^)  are  functions  of  the  two  unknowns  0  and  0.  In 

practice,  sufficient  data  are  available  be  provide  an  overdetermined, 
but  non-linear  system  of  equations.  A  direct,  least  squares  approach  is, 
therefore,  not  possible.  We  resort  to  a  computing  procedure  in  which 
the  unknowns  are  first  approximated  and  then  improved  by  a  series  of 
differential  corrections.  If  second  and  higher  order  terms  are  omitted 
from  the  Taylor  expansion  about  the  point  (0,  0)  which  approximates  the 
solution,  the  equations  of  condition  may  be  written  in  matrix  form  as 

JAW  =  AU.  (28) 


where 


AW  s 


J  s 


du. 


d0  '  B0 


u.  s  cos  X  , 
1  i 


AU  5 


(Von  -  <ui>c 


27 


V 


for  all  values  of  1.  The  subscript  ob  refers  to  observed,  values  while 
the  subscript  c  indicates  computed  values  derived  from  the  approximate 
solution  (6,  0).  Since  the  system  is  overdetermined,  we  obtain  AW  from 
the  least  squares  solution 

AW  -  (J*J)_1  J*  AU  ,  (29) 

where  J*  is  the  transpose  of  the  matrix  J.  Expressions  for  the  evaluation 
of  the  elements  of  J  may  be  found  in  the  Appendix.  Improved  values  for  the 
unknowns  are  obtained  from 

W  +  AW  , 

where 

As  a  matter  of  convenience,  no  subscripts  were  introduced  to  indicate 
iteration.  However,  at  this  point  in  the  computing  procedure,  the 
improved  values  of  W  are  used  as  our  approximate  solution  and  the  process 
is  iterated  until  convergence  is  achieved.  After  convergence,  the  azimuth 
and  elevation  angles,  which  completely  describe  rocket  aspect,  may  be 
computed  from  Equations  (2)  and  (3). 

Frequently,  in  differential  correction  computations,  it  is  difficult 
if  not  impossible  to  initially  approximate  the  unknowns  with  sufficient 
accuracy  to  assure  convergence  to  the  correct  result.  It  is  notable, 
however,  that  the  method  presented  herein  has  been  found  to  converge  for 
any  initial  approximations  of  the  unknown  angles  that  lie  within  the 
permissible  ranges  of  0  to  2tt  radians  for  0,  and  either  0  to  n/2  radians 
or  Tt/2  to  n  radians  for  6.  Tabulated  in  Table  II  are  the  results  of  a 
series  of  computations  for  a  particular  set  of  actual  field  data.  Initial 
approximations  for  0  and  0  have  been  varied  over  the  entire  ranges  of 
these  two  angles.  In  each  instance,  the  computation  converged  to  the 
same  result.  The  number  of  iterations  required  for  convergence  has  also 
been  tabulated.  As  expected,  the  better  approximations  require  fewer 
iterations  and  therefore,  less  computing  time;  but  even  poor  starting 
values  eventually  lead  to  correct  results.  In  practice,  when  aspect  is 
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computed  for  a  sequence  of  points  spaced  at  short  time  intervals,  the 
number  of  iterations,  and  hence  computing  time,  is  held  to  a  minimum  by 
using  the  results  for  previous  points  to  arrive  at  good  initial  approxima¬ 
tions  for  the  unknown  angles. 


TABLE  II 

TEST  OF  COMPUTING  PROCEDURE  FOR  CONVERGENCE 


Computing  Input 

Computing  Output 

Number  of 

a 

€ 

a 

e 

Iterations 

0° 

O 

CT\ 

oc 

18.922° 

86.983° 

5 

90° 

89° 

18.922° 

86.983° 

14 

l8o° 

89° 

18.922° 

86.983° 

11 

ro 

—0 

0 

0 

89° 

18.922° 

86.983° 

9 

0° 

1*5°  ■ 

18.922° 

86.983° 

5 

90° 

U5° 

18.922° 

86.983° 

8 

lp0° 

45° 

18.922° 

86.983° 

8 

270° 

45° 

18.922° 

86.983° 

9 

0° 

1° 

18.922° 

86.983° 

8 

90° 

1° 

18.922° 

86.983° 

13 

l8o° 

1° 

18.922° 

86.983° 

13 

O 

O 
r— 
0 J 

1° 

18.922° 

86.983° 
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DISCUSSION  OF  SYSTEM  ERROR  SOURCES 

All  magnetometer  systems  suffer  from  one  basic  shortcoming;  namely, 
that  elimination  of  magnetic  interference  is  tedious.  However,  the  pay- 
loads  in  which  this  system  has  been  flown  were  designed  primarily  for 
measurement  of  the  geomagnetic  field  and  consequently  the  payload  components 
were  chosen  to  have  small  magnetic  moments.  The  static  magnetic  field 
perturbation  (due  to  permanent  and  induced  magnetism)  at  the  aspect 
magnetometer  sensor  was  measured  for  a  range  of  payload  orientation  and 
found  to  be  less  than  0.1  percent  of  the  earth* s  magnetic  field. 
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Disturbances  due  to  electrical  currents  were  of  similar  magnitude ,  except 
during  the  interval  from  30  to  60  seconds,  when  a  known  disturbance,  due 
to  a  proton  precession  magnetometer,  took  place  during  50  percent  of  the 
measurements.  These  measurements  have  not  been  corrected,  and  consequently 
the  data  scatter  during  this  time  is  large. 

Another  source  of  magnetic  disturbance,  often  neglected  in  the  use 
of  aspect  magnetometers,  is  the  attenuation  and  rotatic  .1  of  the  geomagnetic 
field  within  the  payload  due  to  currents  induced  by  the  field  in  a  spinning, 
conducting  rocket  structure.  The  structure  of  the  nose  cones  under 
discussion  was  fiberglass  and  there  was  no  appreciable  amount  of  metal 
near  the  magnetometers,  so  that  induced  currents  were  negligible. 

Alignment  of  the  magnetic  and  optical  axes  is  easily  handled  by 
standard  machine  shop  practice  to  an  angular  accuracy  greatlv  exceeding 
that  necessary  for  this  application.  The  alignment  of  the  aspect  magne¬ 
tometer  axis  with  the  mechanical  structure  of  the  sensor  is  within  1/4  to 
l/2  degree,  so  that  mechanical  alignment  assures  magnetic  alignment. 

The  systematic  error  in  the  magnetometer  measurement  is  estimated  to 
be  less  than  one  percent,  corresponding  to  an  error  of  l/2  degree  in 
elevation  and  10  degrees  in  azimuth,  or  a  pointing  error  of  less  than  2 
degrees  during  the  vacuum  portion  of  flight.  This  is  the  estimated  maximum 
systematic  error  from  this  source  and  the  actual  value  is  probably 
considerably  less. 

Rocket  aspect  errors  resulting  from  errors  in  the  solar  cell  time 
pulse  have  been  estimated  and  the  results  plotted  in  Figure  7-  A  time 
error  equivalent  to  one  percent  of  the  rocket  spin  period  was  assumed 
and  error  propagation  into  missile  aspect  was  computed  with  data  that 
were  otherwise  free  of  error.  The  results  are  presented  as  a  function 
of  rocket  attitude  with  the  azimuth  angle  plotted  as  the  polar  angle 
and  the  elevation  angle  as  the  radius  vector.  Contour  lines  for  errors 
of  1,  3,  and  5  degrees  have  been  drawn.  The  assumed  time  error  of  one 
percent  of  the  spin  period  is  larger  than  anticipated  for  the  measuring 
system;  yet,  the  resulting  aspect  error  at  higher  elevation  angles  is  not 
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ASPECT  ERROR  RESULTING 
FROM  SOLAR  CELL  TIME  PULSE 


excessive.  It  is  estimated  that  time  errors  in  the  solar  cell  pulse 
will  normally  not  exceed  0.5  percent  of  the  spin  period.  However,  since 
the  aspect  error  varies  almost  linearly  with  the  time  error.  Figure  8 
may  still  be  used  to  estimate  the  effect  of  the  smaller  time  errors  by 
decreasing  the  system  error  in  proportion  to  the  decrease  in  the  time 
error.  The  orientation  of  the  contour  lines  in  Figure  8  is  closely 
related  to  the  direction  of  the  vectors  S  arid  H.  Both  vectors  were 
assumed  to  be  identical  with  those  employed  in  the  reduction  of  the 
measured  data  which  are  presented  in  the  following  section.  It  Is  of 
interest  that  the  minimum  propagation  of  error  occurs  in  the  direction 
of  the  sun  which  was  approximately  30  degrees  west  of  south. 

An  independent  knowledge  of  both  magnitude  and  direction  of  the 
earth’s  magnetic  field  is  necessary  when  only  a  single  aspect  magneto¬ 
meter  Is  employed.  ■  Both  magnitude  and  direction  are  well  known  at  Fort 
Churchill.  Table  III  is  a  list  of  the  measured  and  calculated  values  of 
total  field,  inclination,  and  declination.  Only  a  severe  magnetic  storm 
would  cause  errors  greater  than  those  already  inherent  in  the  system. 

In  the  data  reduction,  an  eccentric  dipole  model  based  on  rocket  measure¬ 
ments^  was  used  to  calculate  the  magnitude  of  the  field.  Declination 
and  inclination  were  considered  constant  throughout  the  rocket  trajectory, 

since  the  variation  of  these  quantities  is  indicated  to  be  less  than 

q 

l/2  degree  from  the  Jensen  and  Cain  model  . 

COMPUTATIONAL  RESULTS 

The  computing  method,  which  has  been  presented,  was  applied  to 

aspect  data  recorded  for  Nike -Apache  Rocket  No.  14.36,  launched  at  Fort 

Churchill,  Canada,  on  7  October  1963-  Reductions  were  obtained  for  most 

of  the  usable  data  and  cover  three  sections  of  the  flight.  The  first 

period,  which  extends  from  5  to  37-5  seconds  after  launch,  is  of  particular 

interest  since  it  is  possible  to  compare  the  results  from  the  aspect 

12 

system  with  missile  attitudes  derived  from  DOVAP  ( DOppler  Velocity 
And  Position)  reductions.  The  aspect  results  for  the  second  period 
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TABLE  III 


TOTAL  SCALAR  FIELD 


Measured  at  Ground  (Proton  Magnetometer) 


Extrapolated  to  Ground  from  Rocket 
Measurement 

.  (9) 

Jensen  and  Cam  Analysis 


INCLINATION 

Measured  at  Ground  (Dip  Circle) 

Jensen  and  Cain 


DECLINATION 

Measured  at  Ground  (Variometer)  (East  of  North) 
Jensen  and  Cain  (East  of  North) 


*10 


5 


Y  =  1  gauss 


61,120  y* 
61,240  y 

61,320  y 

83°  32'  + 
83°  33’ 

2°  35 f  + 
2°  26  1 
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COMPARISON 
AND  ASPEC 
DETERMINA' 


from  265  to  3^2*5  seconds  indicate  that  the  rocket  slowly  nosed  over  as 
it  re-entered  the  effective  atmosphere.  Finally,  after  3 ^6  seconds,  the 
missile  underwent  a  period  of  coning  which  was  continuously  damped  as 
the  missile  gradually  approached  vertical  flight  with  the  nose  downward 
at  375  seconds. 

The  comparison  of  the  aspect  results  with  the  DOVAP  reductions 

yielded  better  agreement  than  anticipated.  The  DOVAP  system  normally 

provides  highly  accurate  determinations  of  rocket  position  and  velocity 

as  a  function  of  time.  However,  missile  spin  introduces  accumulative 

13 

phase  errors  into  the  DOVAP  data  ,  Ordinarily,  these  errors  may  be 
measured  and  the  raw  data  rather  accurately  corrected.  Occasionally, 
as  a  result  of  poor  antenna  patterns,  the  spin  corrections  are  of 
doubtful  quality.  In  the  case  under  consideration,  the  spin  errors  were 
particularly  difficult  to  detect  and  correct.  Therefore,  the  DOVAP 
reductions  very  likely  contain  rather  serious  spin  errors.  These  enter 
the  results  primarily  as  systematic  rather  than  random  errors  since  a 
portion  of  the  spin  correcting  procedure  consists  of  smoothing  data  which 
have  been  corrected  for  the  errors  introduced  by  spin. 

DOVAP  reductions  normally  provide  a  relatively  smooth  trajectory 
which  may  be  differentiated  numerically  to  obtain  rather  accurate  velocity 
components.  If  we  assume  that  the  velocity  vector  lies  along  the  symmetry 
axis  of  the  missile  while  the  latter  is  still  in  the  effective  atmosphere, 
we  may  use  the  DOVAP  velocity  components  to  compute  rocket  attitude  as 
a  function  of  time.  This  has  been  done  for  the  period  from  5  to  37*5 
seconds  and  the  results  have  been  plotted  in  Figures  8  and  9  together 
with  the  reductions  obtained  from  the  aspect  system.  Both  sets  of  data 
show  the  elevation  angle  of  the  rocket  to  initially  decrease  and  later 
increase.  The  latter  event  coincides  with  the  initiation  of  second- 
stage  burning.  While  agreement  between  the  two  systems  appears  to  be 
poor  for  the  azimuth  over  the  early  portion  of  the  trajectory,  it  should 
be  noted  that  the  azimuth  is  not  particularly  significant  when  the 
elevation  angle  approximates  tt/2  radians.  A  more  significant  comparison 
between  the  two  sets  of  data  may  be  obtained  from  Figure  10  which 
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COMPARISON  OF  DOVAP 
AND  ASPECT  SYSTEM 
DETERMINATION  OF 

'*0-  AZIMUTH  ANGLE 


i 


t 
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FIGURE  9 


ANGULAR  DIFFERENCE 


presents  a  plot  of  the  angle  between  the  two  determinations  of  the  rocket 
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maximum  of  about  7  degrees  early  in  flight,  the  two  systems  agree,  on 
the  average,  to  within  1.5  degrees  between  20  and  37.5  seconds.  This  is 
remarkable  since  the  aspect  system  itself  was  designed  for  approximately 
1  degree  accuracy,  and  It  is  estimated  that  the  errors  in  the  DOVAP 
system  are  of  this  order  of  magnitude.  The  extreme  roughness  in  the 
reductions  for  the  aspect  system  in  the  period  from  30  to  37*5  seconds 
was  expected.  It  resulted  from  interference  by  the  polarization  field 
of  a  proton  precession  magnetometer  which  was  also  carried  on  board  the 
rocket  Estimates  of  probable  error  obtained  from  the  least  squares 
solution  for  the  aspect  system  are  on  the  average  about  0.3  degree  in 
elevation  angle  and  6.5  degrees  in  azimuth  over  the  entire  interval  from 
5  to  37*5  seconds.  A  polar  plot  of  the  two  determinations  of  rocket 
aspect  for  the  initial  reduction  period  is  presented  in  Figure  11. 

Azimuth  has  been  plotted  as  the  polar  angle  and  elevation  as  the  radius 
vector.  The  numbered  positions  are  at  five  second  intervals.  In  comparing 
the  reduced  data  for  these  two  systems,  it  should  he  recalled  that  while 
the  DOVAP  results  have  small  random  errors,  they  undoubtedly  contain 
rather  serious  systematic  errors  as  a  result  of  missile  spin. 


The  next  period  of  interest  runs  from  £65  to  3^2.5  seconds.  Plots 
of  both  the  elevation  and  azimuth  angles  are  presented  as  a  function  of 
time  in  Figure  12.  These  graphs  clearly  indicate  that  the  rocket  was 
precessing  at  a  rate  of  approximately  6  degrees  per  second  before  its 
re-entry  into  the  effective  atmosphere  at  approximately  330  seconds. 
Throughout  this  period  of  precessing  motion,  the  elevation  angle  varied 
between  extremes  of  6l  degrees  and  86  degrees  while  the  azimuth  ranged 
from  172  degrees  to  269  degrees.  The  plots  in  Figure  12  show  that  the 
rocket,  following  re -entry;  began  to  nose  over  as  the  elevation  angxe 
decreased  from  85  degrees  to  1*5  degrees  in  the  period  from  330  to  3^2.5 
seconds .  After  a  short  break  of  about  four  seconds  in  the  recorded  data 
of  the  aspect  system,  the  reductions  show  tliat  the  missile  nosed  over 
and  entered  a  stage  of  coning  which  was  gradually  damped  during  the 
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remaining  period  of  observation.  The  aspect  measurements  may  he  easily 
interpreted  from  the  series  of  polar  plots  presented  in  Figures  13A  through 
14a.  As  before,  azimuth  is  plotted  as  the  polar  angle  while  the  negative 
elevation  angle  is  plotted  as  the  radius  vector.  The  terminal  point  of 
each  graph  is  also  the  initial  point  of  the  following  figure.  Flight 
time  has  been  indicated  periodically.  To  aid  in  the  interpretation  of 
the  reduced  data,  continuous  curves  have  been  drawn  through  the  plotted 
points.  With  the  aid  of  these  polar  plots,  curves  showing  the  variation 
of  the  azimuth  and  elevation  angles  with  time  were  derived  and  are 
presented  in  Figures  15  and  16. 


CONCLUSION 

The  aspect  measuring  system,  described  in  this  report,  has  been 
found  to  provide  a  practical  and  inexpensive  approach  for  the  accurate 
determination  of  the  attitude  of  spinning  rockets.  While  the  system 
has  certainly  not  been  precisely  evaluated  for  accuracy,  the  comparison 
of  the  aspect  reductions  with  DOVAP  results  gave  better  agreement  than 
expected.  Recalling  that  the  DOVAP  spin -error  corrections  had  a  low 
confidence  level  for  this  particular  rocket,  it  is  considered  likely 
that  the  errors  in  the  results  for  the  aspect  system  were  smaller  than 
those  in  the  DOVAP  determinations.  On  the  basis  of  these  initial  results, 
the  attitude  sensor  is  estimated  to  have  provided  an  over -all  accuracy 
of  one  to  two  degrees  in  the  determination  of  rocket  aspect.  Moreover, 
it  appears  likely  that  the  sensor  has  a  potential  accuracy  of  better  than 
one  degree. 

The  method  of  solution  has  been  shown  to  be  practical.  It  readily 
lends  itself  to  rapid  computation  on  a  large  scale  digital  computer.  On 
the  average,  one  determination  of  missile  aspect  required  3  to  4  itera¬ 
tions  and  was  computed  at  the  rate  of  0.l6  seconds  per  iteration  on  the 
BRLESC  computer  so  that  a  single  solution  for  rocket  aspect  required 
approximately  0.6  seconds  of  computing  time.  There  are  no  ambiguities 
in  the  results  if  the  direction  of  spin  and  the  sign  of  the  elevation 
angle  are  known.  The  latter  may  be  determined  from  the  magnetometer 
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record  while  the  direction  of  spin  is  usually  a  known  characteristic 
of  the  rocket.  With  ambiguities  resolved,  convergence  to  a  unique  solu- 
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APPENDIX 
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